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Abstract 

The break-up of ^^Be was studied at 41 AMeV using a secondary beam of ^^Be 
from the GANIL facihty on a ^^Ti target by measuring correlations between the 
^"Be core, the emitted neutrons and gamma rays. The nuclear break-up leading to 
the emission of a neutron at large angle in the laboratory frame is identified with the 
"towing mode" through its characteristic n-fragment correlation. The experimental 
spectra are compared with a model where the time dependent Schrodinger equation 
(TDSE) is solved for the neutron initially in the ^^Be. A good agreement is found 
between experiment and theory for the shapes of neutron experimental energies 
and angular distributions. The spectroscopic factor of the 2s orbital is tentatively 
extracted to be 0.46 ± 0.15. The neutron emission from the Ip and Id orbitals is 
also studied. 
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1 Introduction 



The investigation of exotic nuclei has strongly developed during the last two 
decades thanks to the availability of radioactive beams produced by in-fiight 
fragmentation (jH |2l [3]) and structural properties of a large variety of nuclei 
far from stability can be accessed. Recent advances have demonstrated the 
importance of a detailed understanding of nuclear reactions in order to infer 
properties of nuclei such as nuclear shapes, spectroscopic factors or pairing 
correlations (?)• Several reaction mechanisms are used : transfer reactions of 
one or two nucleons (jSl Ej) , Coulomb or nuclear break-up reactions (jH O El |9|) , 
the latter being the focus of this paper. 

With exotic nuclei, we have to deal with low intensity beams and we expect 
an increasing complexity in obtaining nuclear structure information due to 
the reduction of the binding energy. In addition, modification of two-body 
correlations like pairing is expected. 

In this paper, we will demonstrate that the study of nuclear break-up around 
the Fermi energy, conjointly with the development of accurate nuclear reaction 
models ( 1T0| [TT| |12j [13]), can provide an alternative path to get information on 
structural properties in stable and weakly bound exotic nuclei. In this respect, 
^^Be appears as a good case for study as it is a one-neutron halo nucleus and 
is abundantly produced. Its ground state is known to be a mixture of a cold 
^''Be core coupled to a halo neutron in a 2si/2 state with an excited ^°Be in 
a 2+ state coupled to a ld5/2 neutron leading to a 1/2 ground state with a 
positive parity. 

\GS > = a|0+®2si/2> + ^|2+®ld5/2> 

The determination of the relative importance of these two states was addressed 
through several experimental and theoretical approaches and is related to the 
spectroscopic factors S2s=a^ and Sirf=/?^. ^^Be was extensively studied during 
the last decade as one of the first observed halo nuclei (jl]) . It has revealed many 
of its properties through several different types of reactions such as break-up 
experiments ([71 [TH [151 [1B|) or transfer reactions ([5i). Coulomb dissociation has 
been shown to play an important role when ^^Be interacts with a heavy target 
such as Au, as the halo neutron is bound by only 500 keV and easily separates 
from the core. 

For lighter targets such as Be, the nuclear break-up is expected to take over 
(ITTll . It has been shown in a time dependent non-perturbative calculation 
(TDSE) ( ITOl [TTjl that nuclear break-up presents a typical emission pattern, 
where the neutron and the remnant ^°Be core are expected to be emitted in 
opposite sides of the beam direction. This property has already been observed 
with stable nuclei and was called "towing mode" (1T5|) . The pattern of this 
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mechanism, leading to the emission of the particle at large angle, is expected 
to be sensitive to the initial wave function such as its angular momentum and 
could be used to shed some light on the structure of the nucleus. 

In this paper, we will report on an experiment performed at the GANIL facil- 
ity on the break-up of ^^Be where a ^°Be fragment was detected in coincidence 
with a neutron and gamma rays. After a description of the experimental ap- 
paratus in section [2] and of the theoretical framework in section [3| we will 



present the neutron angular distribution in paragraph 4.1 and demonstrate 
the characteristic anti-correlation between the neutron and the ^°Be in para- 



graph 4.2 The gamma energy distribution will be presented in paragraph 4.3 



Paragraphs |4.4 and 4^ will be devoted to the neutron energy spectra in two- 



fold and three-fold events respectively. By comparing the experimental data 
with the time dependent Schrodinger calculation described in section |3| the 
possible extraction of spectroscopic factors of the ground state of ^^Be will be 
discussed in section |5l 



2 Experimental set-up 



The reaction ^STi(i^Be,^°Be+n) at 41 MeV per nucleon was studied. The ^^Be 
beam was produced by fragmentation of a 10 /iA i3q6+ pj-imarj beam onto 
a 6128 fim thick carbon target placed in the Superconducting Intense Source 
for Secondary Ions device (SISSI) of the GANIL facility. A 600 /im thick 
aluminum degrader was set in the alpha spectrometer and a pure ^^Be beam 
was transported to the SPEG vault with a small angular divergence of 0.2". 
The resulting beam intensity of 80000 pps was monitored by a drift chamber 
placed before the target. The beam size on the target was ±6.4 mm giving 
an emittance after the break-up of more than 400 vr.mm.mrad and a loss in 
detection efficiency of the SPEG (|T9|) spectrometer, which has a 70 7r.mm.mrad 
nominal acceptance, of about a factor of 5. However, with the small angular 
divergence, no beam tracking was needed to extract the scattering angle of 
^"Be and hence to perform angular correlations between the ejectile and the 
neutrons. After breaking up on a 113 mg/cm^ Ti target, the ^°Be core was 
analyzed by the energy loss spectrometer SPEG positioned at degree and 
covering ±2° in both horizontal and vertical directions. 74 BaF2 scintillators 
covered 84% of the total solid angle around the target to detect the gamma 
rays and 23 neutron detectors were utilized and positioned as shown in the 
experimental set-up presented in Figjl] 

Two drift chambers, an ionization chamber and a plastic detector forming the 
standard detection in the focal plane of SPEG were used for identifying the 
fragments and infer their momenta and scattering angles. An identification 
plot of the fragments produced after the Ti target is constructed by plotting 
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Fig. 1. Experimental set-up : The 74 gamma rays detectors were placed below 
and above the horizontal plane (37 detectors are shown on the plot). 23 neutron 
detectors were positioned around the target and the ^'^Be fragments were detected 
after analysis by the spectrometer SPEC 




Fig. 2. Identification of the fragments produced in the Ti target and after the SPEG 
spectrometer. 

the energy loss in the ionization chamber versus the arrival time of the ion in 
the plastic detector (proportional to the M/Q of the particle since the Bp is 
fixed by the spectrometer) and shown on Fig|2} The spectrometer field was 
set to obtain the ^°Be in the middle of the focal plane. 
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Fig. 3. Typical time of flight spectrum for the neutron detectors. 

The neutrons were detected by 11 detectors of the NordBall array (2D) plus 
8 cylindrical detectors of either 2 or 6 inches thick from the Uppsala group 
and 4, 30 cm long cylindrical detectors from Orsay. These neutron detectors 
of BC501 type were placed around the target at a distance varying from 1.4 
meters to 3.5 meters (see Figjl]) and at angles running from 30° to 95" in the 
laboratory frame. The identification was given by pulse shape discrimination 
by integrating the signal in two 400 ns wide gates, the first starting at the 
beginning of the signal and the second shifted by 30 ns. Energies were obtained 
from the time of flight (ToF) measured between the RF of the beam and the 
neutron arrival time in the detector. A typical ToF spectrum is presented in 
FigjS} 



74 BaF2 detectors of the Chateau de Cristal ensemble were used to detect 
the gamma rays. Detector efficiency and energy calibrations where performed 
with ^°Co and ^'^^Cs radioactive sources. GEANT3 was used to simulate the full 
detector and gave an overestimation of the experimental efficiencies obtained 
at 661 keV (^^^Cs), 1173 and 1332 keV C^^Co) by 23%, as already observed in 
ref . (12T|) . The detection efficiency of gamma rays of higher energies was then 
obtained by the GEANT3 simulation code (release 3.21). Furthermore, as 
the cells were generally not directed toward the target and were overlapping 
in angle, GEANT3 was also used to infer the effective angle at which each 
detector measured the gamma rays coming from the target position. These 
angles, that were different from the average geometrical angle of each detector, 
were necessary for Doppler correction of the energy of the gamma rays emitted 
by the ejectile. As will be shown in paragraph 4.3, three peaks corresponding 
to gamma rays from ^°Be were clearly observed. Whenever these peaks could 
be seen for a single detector in coincidence with a ^°Be, their positions were 
used for the energy calibration at high gamma ray energies. The energy range 
of interest in our case is about 3.37 MeV in the ^"Be frame as it corresponds 
to the decay of the 2+ state to the OS (see paragraph 4.3) which corresponds 
to 4.5 MeV in the laboratory frame at very forward angles. 



5 



3 TDSE Calculations 



During the past decade several approaches have been developed to describe 
the Coulomb and nuclear breakup at and above the Fermi energy. Among 
them, the Eikonal approximation, which takes advantage of the high velocity 
of the projectile, has encountered a large success in particular to interpret 
recent knock-out reactions (see for instance the recent review (jB])). It is how- 
ever worth mentioning that the "towing- mode" which is part of the diffractive 
dissociation channel mechanism is very specific. Indeed, as noted in ref. (ITTj) . 
description and understanding of the "towing mode" goes beyond calcula- 
tions within the Eikonal approximation. On opposite, models based on the 
resolution of TDSE equations ([TnilIIl[l3l[l2l|22l[23l|2ll[25l|26|) properly 
account for such a continuum emission. The great advantage of TDSE tech- 
niques is to treat non-perturbatively the nucleon emission to the continuum 
and the interferences between Coulomb and nuclear break-up. On the other 
hand, TDSE method often neglect interferences due to the mixing of different 
impact parameters that are generally treated properly in the Eikonal model. 
Only recently, methods have been developed that combines the advantages of 
both Eikonal and TDSE framework (JQj 



Here, we restrict ourselves to the simplified situation where a time dependent 
Schrodinger equation (TDSE) has been solved for a one-nucleon wave function 
initially in a ^°Be core and a comparison with the experimental data will be 
presented in section [5] to extract spectroscopic factors. In the present section 
we will briefly describe the calculation and more details can be obtained from 
ref. ([101 HID. 

The initial wave function of the halo neutron was deduced from the solution 
of the single-particle Hamiltonian that writes 

= + Vs.(f) (1) 
where the ^"^Be potential was taken as a Wood-Saxon 

Vseir) = %^ (2) 

1 + e "0 

where the diffuseness = 0.75 fm and the radius Rbs = 1.27.11^/'^ =2.82 fm. 
Numerically the initial wave-packet is obtained by diagonalizing the one-body 
Hamiltonian in spherical coordinates in a sphere of 30 fm radius with a space- 
step of Ar = 0.02 fm. The depth Vq was taken to be -46.025 MeV to finally 
obtain a 2s wave function bound by 0.503 MeV. 

The wave function calculated in spherical coordinates for ^^Be is then trans- 
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formed in Cartesian coordinates for the dynamical evolution and an appro- 
priate interpolation is made. Special attention has been given to the purity 
of the ground state ^^Be neutron halo by performing an additional imaginary 
time evolution on the Cartesian network, after which the 2s state was bound 
by -0.503 MeV. 

The dynamical evolution of the system is given by the following single particle 
Schrodinger equation; 



where Vr and Vbc are the target and the projectile time- dependent potentials. 
The positions rxit) and rBe{t) correspond to the target and the projectile 
respectively. 

The nuclear potential of the Ti target is taken to be a Wood-Saxon potential 
Vrir) with a diffuseness a equal to 0.75 fm, a radius R = 1.17^4^/^ fm (A 
being the mass number of Ti) and a depth Vq of -41.2 MeV adjusted to obtain 
the experimental binding energy of the last neutron in a If state which is 
-11.6 MeV. The sensitivity to different parameters of the calculation will be 
discussed in the following. 

In order to take into account the ^"^Be core displacement, the evolution is 
performed for each time step using the classical Coulomb trajectory for the 
center of mass motions from a distance of -400 fm along the initial velocity 
axis between the projectile and the target and equal to the impact parameter 
b on the perpendicular axis. The Runge Kutta numerical method is applied 
for the trajectory calculation. The spht operator method (1281) is used to solve 
the time dependent Schrodinger equation with a time step of 1.7 fm/c on a 
mesh of 64x64x64 (fm^) with a step size of 0.5 fm. This method precludes the 
inclusion of a spin-orbit term. 

The calculation is then performed for several impact parameters with either 
a sharp cut-off or a weight that corresponds to a survival factor mimicking a 
probability that the core remains intact during the collision (:29). This prob- 
ability writes: 



where = lA{Al^^ + (Ap — lY^^) and = 0.6 fm. For the considered reaction, 
we have = 8.1 fm. The maximum impact parameter, corresponding to the 
impact parameter where the fraction of the wave function emitted saturates, is 
110 fm in the case of a Ti target. At this impact parameter the fraction of wave 
function emitted is around 0.04% and only the noise inherent to the numerical 




+ Vrif - TTit)) + Vseif - rBe{t))}\(p) 



(3) 




(4) 
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Fig. 4. Probability density of the neutron wave function after the scattering of 
^^Be onto a Ti target at an impact parameter of 8 fm. The plot is presented in 
logarithmic scale and summed over the coordinate perpendicular to the reaction 
plan. Also presented is the trajectory of the ^^Be and the initial position of the 
target. 



method remains. No larger impact parameter is included in the cross section 
estimate. 

The result of the calculation is presented in Figj4]as a probability density plot 
of the neutron for a collision with a Ti target at an impact parameter of 8 
fm. We observe that a fraction of the wave function is emitted at large angle. 
This is what is called the "towing mode", already experimentally observed 
with stable beam (IT^ . 

Besides the particle emitted to the continuum, the rest of the wave-function 
either remains in the projectile potential or is transfered to the target. Since 
in that case, the particle is not emitted, these contributions are removed to 
compare with the break-up channel measured in the experiment. In practice, 
this is performed by adding an imaginary part to the target and projectile 
potentials at the end of the calculation. 
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Fig. 5. Angular distributions for neutrons of kinetic energy above 26 MeV. Open 
dots are data from ref. (fTi|l and black dots are our data. 

4 Experimental results 



4-1 Neutron angular distribution in coincidence with ^^Be 



FigjSj shows the neutron angular distribution in coincidence with ^''Be ejectiles 
and compared to the experimental data of ref. (11311 with the same neutron 
energy threshold of 26 MeV. The data are corrected for the neutron and the 
spectrometer detection efficiencies. At neutron angles above 30°, due to the 
use of a spectrometer that allowed for the removal of the incident beam before 
the fragment was detected, much higher statistics was gathered allowing for a 
more precise analysis as will be shown in the following. 



4-2 Angular correlation between the ^^Be and the neutron 



A previous analysis (1T7I) of the data of ref.(jllj) has shown that the large angle 
neutron emission is mainly due to small impact parameters (b < 15 fm) for 
which the halo wave function is diffracted by the nuclear potential of the 
target. In such a case, the TDSE calculation predicts that the ^''Be and the 
neutron are preferentially emitted on opposite sides of the beam, giving rise to 
a specific angular correlation as shown in Fig|6j It shows that when the core is 
scattered on the right side of the target, most of the neutron wave function is 
present on the opposite side. In the experiment of ref. (JEj), the angle of ^''Be 
were not measured, precluding the investigation of such correlations. 

In this experiment, such a correlation could be extracted since the scattering 
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Fig. 6. Neutron angular distribution (after efficiency corrections) for neutron en- 
ergies above 26 MeV measured in coincidence with ^'^Be ejectiles. Black dots are 
for events where the neutron and the ejectile were measured on opposite sides of 
the beam (|$r-e«| > 90°) whereas the open circles are the events where the neutron 
and the fragment have been detected on the same side of the beam (|$re/| < 90°). 
Corresponding TDSE calculation performed for a 2s wave function are represented 
as the plain and the dashed lines, respectively. 

angle of the ejectile was measured in the SPEG detection system. The exper- 
imental neutron angular distributions are shown after efficiency correction in 
Fig|6]for events where the ejectile and the neutron were detected on opposite 
sides of the beam axis (black circles, $re/ > 90°) and for events where both the 
ejectile and the neutron were detected on the same side of the beam (empty 
circles, ^rei < 90°). We observe an enhancement of about a factor 4 for the 
neutron cross section when detected in coincidence with a ^°Be on the opposite 
side of the beam with respect to the same side of the beam. The solid curve 
results from the calculation for a neutron and a core on opposite sides of the 
initial ^^Be direction arbitrarily normalized to the data points (black circles). 
The dashed line corresponds to a core and a neutron flying on the same side. 
Presented calculations were performed for a 2s wave function but other angu- 
lar momenta show no difference in neither the neutron angular distributions 
nor the relative asymmetry. The relative normalization of these two curves is 
given directly by the calculation. In the region between 30° and 50°, where the 
cross section is large, the calculation reproduces very well the experimental 
asymmetry. At larger angles the actual asymmetry is slightly smaller than the 
prediction. Overall the anti-correlation between ^°Be and neutrons, which is a 
very detailed feature of the reaction mechanism, is remarkably well accounted 
for. 

We thus conclude that the break-up mechanism is well understood and de- 
scribed by our dynamical calculation and that the emission of neutrons at 
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large angles is driven by the nuclear interaction of the neutron with the target 
through the so-called "towing mode" mechanism. 



4.3 Gamma spectrum 



In order to access the final states of the ^°Be core after the break-up, coin- 
cidences were measured with the gamma rays emitted during the reaction. 
Figure [7] presents the gamma ray spectrum, to which a Doppler correction 
was applied to account for the velocity of the ejectile, recorded in coincidence 
with a ^°Be. It shows three peaks located around 2.6, 2.9 and 3.4 MeV, which 
correspond to the known transitions of ^°Be shown in the level scheme dis- 
played in the right panel. To estimate the background, a run with a ^'^O beam 
was performed and the gamma spectrum extracted in coincidence with a ^^N 
detected in the spectrometer. In this reaction, no gamma rays above 1 MeV 
are expected. The whole Chateau de Cristal configuration was described in the 
GEANT3 code to infer its response to these three observed peaks as well as the 
direct decay from the 1~ state to the GS which gives rise to a gamma ray of 
5.96 MeV energy (130|1 . A fit was performed including these four components 
plus the background, with free normalizing factor for each, and reproduces 
well the experimental spectrum above 1.5 MeV. 

The level scheme presented in FigjT] shows that the gamma ray of the decay 
from the 2+ state to the GS is part of the cascade of two other decays from 
the 1~ and the 2~ states. Therefore, the observation of a 3.37 MeV gamma 
ray signifies that either we started from a 2"*" state with a neutron hence in a 
Id state, or we are measuring the decay from the 1~ or the 2~ states through 
the 2"*" state and the emitted neutron stemmed from the ^"Be core, hence from 
a lp3/2 or a lsi/2 state. Therefore, the number of 2+ to 0^ gamma rays should 
be equal to the number of cascades from the 2~ and the 1~ states plus the 
number of events feeding directly the 2"*" state. 

From the fit to the data presented in FigjT] we extracted that 2/3 of the 
3.37 MeV gamma ray were the result of a cascade from the 2~ or the 1~ 
state, whereas 1/3 of the total number of 3.37 gamma rays were due to a 



direct feeding of the 2"^ state. This ratio will be used in paragraph |4.5| when 
analyzing three-fold events. 



4-4 Neutron kinetic energy spectra - 2 fold events 



The neutron energy spectrum for two-fold events where a neutron was de- 
tected above 30° in the laboratory frame in coincidence with a ^°Be (black 
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Fig. 7. Left: Gamma energy spectrum in coincidence with a ^'^Be and fitted by the 
contributions of 4 peaks with the GEANT3 code. Right: Level scheme of ^^Be. 
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Fig. 8. Neutron energy spectrum in coincidence with a ^'^Be and for angles above 
30" (black dots). Open circles are the data in coincidence with a ^^Be. 



dots) is presented in FigjSj It shows a peak centered at about 40 MeV and a 
contribution at low energy. 



The angular distribution for neutrons above 26 MeV has already been shown 
in FigjS] whereas the angular distribution for neutrons below 15 MeV is pre- 
sented in FigjOj The latter shows an isotropic behavior in the laboratory frame 
indicating that these neutrons are not coming from the moving projectile but 
rather from the target at rest. A way of estimating the shape of the energy 
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Fig. 9. Angular distribution of neutrons below 15 MeV of kinetic energy. 
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Fig. 10. Neutron energy spectrum in coincidence with a ^"^Be and subtracted for the 
neutrons emitted by the target. 



spectrum of the neutrons emitted by the target is to extract the neutron spec- 
trum in coincidence with a ^^Be detected in the focal plane. Here, apart from 
a pick-up break-up component which should be very small, one ensures that 
the detected neutron arises from inelastic excitation followed by decay of the 
target. This spectrum is presented in Figjs] (open circles) and exhibits only 
a low energy neutron component. We now subtract this component from the 



spectrum of FigjS] (black dots) and the result is presented in Fig, 10 The high 
energy peak now dominates the spectrum with a small tail extending to lower 
energies. 
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Fig. 11. Neutron energy spectra in coincidence with a ^^Be and a gamma ray of 
energy between 1 and 4 MeV. Left: the black dots are the total exclusive data and 
the open circles are the background coming from the target emission and normalized 
to best fit the black dot spectrum. Right: result from the subtraction of the target 
emission. 

4-5 Three-fold events 



We present on Fig, 13 the neutron energy spectrum in coincidence with gamma 



rays between 1 and 4 MeV. We estimate from the GEANT3 simulation that the 
efficiency for detecting the decay of the 2+ state is 50% in this energy range. 



Furthermore, according to paragraph 4.3, 1/3 of the 2"*" to the O"*" gamma-rays 
come from direct feeding of the 2+ state. Besides the low energy component, 
a broad structure is observed, located around 30 MeV. In the same ways as 



in paragraph 4.4, the subtraction of the neutron component coming from the 
target (open circle) is then performed and leads to the right panel of Fig, 11 
where only the 30 MeV structure remains. 

In order now to obtain a neutron energy spectrum for events where no gamma 



ray was emitted we subtract the three-fold spectrum (Fig, 11 right) from the 



inclusive spectrum of Fig, 10 The result is presented in Fig, 12 as black dots. 
For this last subtraction, no normalization was performed as the detection 
probabilities were taken into account. A peak clearly shows centered about 40 
MeV. 



Comparison with the TDSE calculation and possible extraction 
of spectroscopic factors 



The TDSE calculation was performed for several initial neutron wave functions 
and results are shown in Fig, 12 For the Id wave function, the ^°Be potential 
is no longer spherical. Indeed, the 2+ state of ^°Be corresponds to a deformed 
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Fig. 12. Neutron energy spectrum when no gamma ray is emitted. Curves are the 
TDSE calculations for a neutron initially in a 2s, Ip and Id state. 



state and was found (jST]), through inelastic scattering of protons, to have 
a deformation parameter /32 = 0.74. To be able to treat the emission from 
a deformed core, the TDSE model has been adapted by slowly deforming 
the core potential during the imaginary time procedure. The iterations are 
stopped once the deformation corresponds to the P2 value of the core and the 
imaginary time procedure has converged. The depth of the obtained deformed 
potential was adjusted to obtain a binding energy of 3.8 MeV, the sum of the 
^^Be binding energy and the 2"*" energy of ^"Be. We performed then the time 
dependent evolution of the calculation with three orientations of the main 
axis, along X, Y and Z to take into account the random orientation of the 
nucleus when the collision occurs. The sum of the calculations for the Id state 



is shown in FigjT2j TDSE calculations were also performed for the a Ip initial 
wave functions with a binding energy of 6 MeV and the result is shown in 



Fig- 12 The calculation also predicts that the contribution from a Is initial 
wave function is negligible compared to the one of the Ip and is not included 
in the following analysis. 



5. 1 Contribution of the 2s orbital 



The experimental neutron energy spectrum for which no gamma ray has been 
emitted is fairly well reproduced by the TDSE calculation performed with 
the 2s wave function. We clearly see in Fig, 12, that the TDSE result is highly 
sensitive to the initial angular momentum of the wave-function. It is also worth 
mentioning that the binding energy of the emitting level influences the kinetic 
energy spectrum. For instance, if the 2s wave-function would have been bound 
by 2.75 MeV, the distribution would have been peaked around 36 MeV instead 
of 40 MeV. 
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Extraction of spectroscopic factors is performed by comparing experimental 
absolute cross-sections with calculated ones. Here, arbitrary normalizations 
were used for the calculations except for the 2s state for which a best fit 
to the data was performed. The calculation, assuming that a 100% of the 
neutrons are coming from the 2s state, provides, with the minimization 
method to best reproduce the data, directly the spectroscopic factor for this 
state. As we have mentionned in section |3j the choice of the minimal impact 
parameter can influence the result. If a sharp cut-off at hmin = 8.5 fm is used, 
the deduced spectroscopic factor is 0.46 while the prescription of ref. (129|1 
leading to a smooth cut-off gives a spectroscopic factor of 0.42. 



5.2 Contribution of the Ip and Id orhitals 



The same comparison was performed for the three-fold events where a gamma 



ray is emitted. Experimental data are shown in Fig, 13, In that case it is 



clear that the data can not be reproduced by the calculated 2s component 
which peaks at 40 MeV whereas the data are centered at about 30 MeV. 
This precludes the two-step process in which a 2s neutron would be emitted 
together with an excitation of the core in a 2"*" state as proposed in ref. fj32jl . 
Both the Ip and Id contributions present a shape similar to the data. However 
a constraint is that we expect to measure twice as many neutrons coming from 
the Ip state than from the Id state. Therefore we fixed this ratio between the 
two calculations and normalized the weighted sum to best reproduce the data. 
This is presented in FigjT3] and the respective spectroscopic factors obtained 
are : Sip = 3.9 ± 1.5 and Sid = 0.50 ± 0.20 where the quoted errors are purely 
statistical. Here, for the Id contribution, we have supposed a sharp cut-off 
between 7 fm and 9 fm depending on the relative orientation of the nuclei. 



5.3 Sensitivity to input parameters on extracted spectroscopic factors for the 
2s orbital 



While the shape of the distribution is rather insensitive to model parame- 
ters, the absolute cross sections appear to be rather sensitive to the inputs of 
the calculation. A preliminary study of the influence of the parameters can be 
found in ref. (1301) . Similarly to knock-out (jSD or transfer ( 13 3 p reactions, param- 
eters should be carefully optimized to be able to infer spectroscopic factors. 



In Fig, 14, we further illustrate the sensitivity of the extracted spectroscopic 



factor by varying the impact parameter cut-off and target potential depth. 

Depending on whether we assume a sharp cut-off in the impact parameter 
at bmin = 8, 8.5 and 9 fm or a weighted sum as described in section [3] the 
spectroscopic factor extracted from the calculation runs from 0.39 to 0.54 for 
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Fig. 13. Neutron energy spectrum in coincidence with a gamma ray of energy be- 
tween 1 and 4 MeV. Curves are the TDSE calculations for a neutron initially in a 
Ip and Id state. The sum of the two contributions is also presented. 
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Fig. 14. Variation of the spectroscopic factor with the minimum impact parameter 
and the absolute value of the perturbating nuclear potential depth. The ellipse 
indicates values corresponding to the ^^Ti potential depth leading to a 1/ bound 
neutron at the proper energy. 

the 2s state when other parameters remain unchanged. We also display in 



FigjT4]the effect of a change in the target potential depth for the sharp cut-off 
bmin = 9 fm and the weighted case. The deeper the potential, the lower is the 
extracted value of the spectroscopic factor. However, it should be noted that 
the prescribed value of 41.2 MeV is optimal to properly describe the target 
binding energy and can not be considered as a free parameter. 

The sensitivity of the calculation has been further tested with respect to the 
radius of the Wood-Saxon potential of the ^^Be nucleus. We used ro=1.20 
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instead of 1.27 and reajusted the depth of the potential to obtain a 2s state 
bound by 0.503 MeV. The extracted spectroscopic factor for bmj„=8.5 fm 
turned out to be 0.47 instead of 0.46. 

Finally, we took the nuclear potential used in ref. (ll7p that has a surface term 
according to N.Vinh Man's prescription ( 13^ . In that case we obtained a lower 
value, 0.42, probably due to the fact that the wave function is on the average 
positionned more out-side the nucleus core and therefore a larger fraction of 
it is towed by the nuclear potential. 

It is important to note that the present TDSE calculations are made with 
several approximations. First, a simple Coulomb trajectory for the target and 
projectile is supposed. Since we are considering here rather small impact pa- 
rameters, a proper treatment of the nuclear plus Coulomb trajectory is a 
priori required. However, to see the influence of the trajectory on the result, 
a straight line trajectory was used and the extracted spectroscopic factor was 
then 0.48 instead of 0.46. It can then be infered that taking into account a 
more realistic trajectory, with the Coulomb and the nuclear potentials, would 
affect only slightly the spectroscopic factors. The second simple approximation 
made regards the fraction of the neutron wave-function initially in the pro- 
jectile and transferred to the target. This transfer is actually poorly treated 
in the current TDSE calculation for two reasons. First, the effect of the Pauli 
principle is neglected. Second, a proper treatment of this mechanism would 
also require to add an imaginary potential to the Ti potential. In our model, 
such an imaginary potential is only added a posteriori to remove the fraction 
of the transferred wave-function when comparing with experimental data. 



6 Discussion 



The value of the spectroscopic factor obtained in this paper is plotted 



in Fig, 15 (adapted from ref. (jTj)), and compared with different experimental 
data and shell model calculations. We see that our result is below the average 
value of other observations but it should be noted that our error bar is rather 
large. During the analysis we had to correct for the small acceptance of the 
spectrometer (see section |2|. Due to the shape of the cut-off, this correction 
brings an error which we can roughly estimate to be between 10 and 20%. 
This is based on FigjS] where a comparison is made between the present data 
and a previous measurement. A rough estimate of the errors on the spec- 
troscopic factor accounting for the acceptance uncertainty, the experimental 
statistical error and the uncertainty on cut-off in impact parameters leads to 
a spectroscopic factor equal to 0.46 ± 0.15. 

As discussed in paragraph |5]3 TDSE calculations performed to extract spec- 



18 



to.a 1 



0,6 - 



o 

o 

u 
« 

I 0,4 
« 

g 

1 0,3 
CO 



Transfert 



ii 



Break-up reactions 



i 



t 



Present 
work ■ 



SM[421 
SM [41J 



Fig. 15. Spectroscopic factors of the |0+ 2si/2) deduced from transfert reactions 
daaiMHSTI ID and nuclear and/or coulomb break-up (1381 ITSj l39l liOt [71 [T6|) . The two 
lines indicates the prediction of two different shell model (SM) calculations (|4H H2 |) . 



troscopic factors are mostly sensitive to the cut-off at low impact parameters. 
In order to obtain a spectroscopic factor of 0.8 for the 2s state we shall use 
a sharp cut-off of about 10 fm which seems quite high considering the two 
nuclei. 

The values extracted in the transfer reaction by Fortier et al. (|5j) (0.84), Wind- 
field et al. (l371) with the same data (from 0.67 to 0.8 depending on the optical 
potential used) and from break-up reactions at relativistic energies of ref. 
(0.61(5) or 0.77(12) depending on whether Coulomb or nuclear break-up is 
selected), show the influence of the experimental method and of the model 
used to describe the collision. Furthermore a value as low as 0.35 has also 
been predicted for transfer reactions (1361) when including the recoil excitation 
and the breakup of ^^Be and the variational shell model gives a value of 0.55 



for the spectroscopic factor of the 2s state (l^fsee dotted line of Fig, 15) 



As for the Id state, the spectroscopic factor found of 0.50±0.2 is much higher 
than the value obtained in the transfer reaction of ref.(E]) but again close to 
the value of 0.40 given by the variational shell model (IH]) that insists on 
the indispensable coupling between the two components of the ground state 
of ^^Be to describe correctly the parity inversion of the l/2~ and the 1/2"*" 
states. 



7 Conclusion 



In this paper we have reported on an experiment studying the break-up of ^^Be 
performed at the GANIL facility. The measurement of the ^°Be remnant in 
coincidence with neutrons emitted at large angles (> 30°) has allowed to select 
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the nuclear break-up ("towing mode"), and has led, with a comparison to a 
TDSE calculation, to the extraction of spectroscopic factors for the ground 
state of ^^Be. 

A value of 823 — 0.46 ± 0.15 is extracted for which a large part of the un- 
certainty comes from the sensitivity to models parameters and experimen- 
tal errors (statistics + systematics). We also obtained Sip = 3.9 ±1.5 and 
Sid = 0.50 ± 0.20 where the error bars only come from the experimental 
statistics. 

Cross sections for the nuclear break-up are of the order of several tens of 
millibarns in the case of weakly bound nuclei such as halo nuclei. If a good 
knowledge of the minimum impact parameter could be achieved and as clean as 
possible experimental data are obtained, we believe that this mechanism could 
be competitive for the study of nuclear spectroscopy of very rare isotopes. 

A study of borromean nuclei such as ^He and ^^Be through nuclear break-up is 
currently in progress. In these cases the measurement of both the neutrons of 
the halo is performed and compared to an appropriate description of the two- 
neutron wave function with a TDSE calculation, extended in such a way to 
include the correlation between the neutrons. Information on the two-neutron 
configuration could then be extracted. 
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